Some topics concerning the theory of singular 

dynamical systems 

Narciso Roman-Roy* 

Departamento de Matematica Aplicada IV. 
Ed. C-3, Campus Norte UPC. C/ Jordi Girona 1. 08034 Barcelona. Spain 



Abstract 

Some subjects related to the geometric theory of singular dynami- 
cal systems are reviewed in this paper. Li particular, the following two 
matters are considered: the theory of canonical transformations for 
presymplectic Hamiltonian systems, and the Lagrangian and Hamil- 
tonian constraint algorithms and the time-evolution operator. 
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1 Introduction 



The aim of this paper is to carry out a brief review of several topics concern- 
ing the theory of autonomous singular dynamical systems, from a geometri- 
cal perspective. In particular, our interest will be focused on two subjects, 
namely: the theory of canonical transformations for singular systems, and the 
problem of the compatibility of the dynamical equations of Lagrangian and 
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Hamiltonian singular systems; more precisely, the analysis of the Lagrangian 
and Hamiltonian constraint algorithms and their relation. 

The article is based on the developments made on the references ^Hj and 
for the theory of canonical transformations, and jl], ^21, ^Hj, ^H] 
and [2ni for aspects related with the constraint algorithms. Thus, we will 
refer to these articles for more details on all these results. 

All manifolds are real and C°°. All maps are C°°. Sum over crossed 
repeated indices is understood. Throughout this paper i{X)u denotes the 
contraction between the vector field X and the differential form u, and L(X)u; 
the Lie derivative of u with respect to the vector field X. 

2 Canonical transformations 

2.1 Presymplectic Locally Hamiltonian Systems 

Definition 1 A presymplectic locally Hamiltonian system (p.l.h.s.) is a 

triad {M,uj,a), where (M, tu) is a presymplectic manifold, and a G Z^{M) 
(i.e., it is a closed 1-form), which is called a Hamiltonian form. 

A p.l.h.s has associated the following equation 



If X exists, it is called a presymplectic locally Hamiltonian vector field as- 
sociated to a. Nevertheless, in the best cases, this equation has consistent 
solutions only in a submanifold jc'- C •—>■ M, where there exist X G X(M), 
tangent to C, such that 



Furthermore, the solution X is not unique, in general, and this non-uniqueness 
is known as gauge freedom. In general (C, lo^ = JcUj) is a presymplectic man- 
ifold which is called the final constraint submanifold (f.c.s.). 

The following theorem gives the local structure of p.l.h.s. (see 0): 
Theorem 1 Let {M,uj,a) be a p.l.h.s., and jc'- C •—>■ M the f.c.s. Then: 



i{X)iJ = a 



X G X(M) 



[t{X)u; - a]\c = 



(1) 



N. Roman-Roy, Some topics concerning the theory of singular systems 3 



1. There are a symplectic manifold (P, fi) and a coisotropic embedding 
ic '■ C "-^ P such that uj^ = t^Q. 

2. For every X G X{M), tangent to C, solution to (OJj, there exists a 
family of vector fields X{P,C) C X(P) such that, for every G X{P,C), 
(a) X^ are tangent to C, (b) X^\c = X\c, and (c) X^ are solutions to 
the equations i{X^)Q = ap + ^, where ap G Z^{P) satisfies i^ap = j'^a, 
and C, G Z^{P) is any first-class constraint form (i.e., t^C, = 0, and the 
Hamiltonian vector field associated with ^, X^ G X(P), is tangent to C). 

3. The coisotropic embedding ic and the family X(P, C) are unique, up to 
an equivalence relation of local symplectomorphisms reducing to the identity 
on C . 

{P,C,Q) is the coisotropic canonical system associated to {M,u!,a). 

2.2 Canonical Transformations for p.l.h.s. 

Let Xih{C) be the set of locally Hamiltonian vector fields in {C,ujc)', that is, 
XihiC) = {Xc G X(C) I UXc)uJc = 0}. 

Definition 2 Let {Mi,uji,ai) (i = 1,2) be a p.l.h.s., with f.c.s. : Ci — >■ 

Mi, and uo^ = j* uoi, such that dim Mi = dim M2, dim Ci = dim C2, and 
rankcijp^ = rankcj^^. A canonical transformation between these systems is a 
pair with $ G Diff (Mi, M2) and G Diff (Ci, C2), such that: 

1. $ o = j^^ o (j); that is, we have the commutative diagram 

Ml M2 

A — ^ ft 

2. 4>.(Xih{Ci)) c X,i(C2). 

The generalization of Lee Hwa Chung's theorem to presymplectic mani- 
folds allows us to prove that (see [IHl): 

Proposition 1 Condition 2 is equivalent to saying that there exists c G M 
such that jI;^{^*uj2 - cui) = (j)*u^ - cujc =0- 
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c is called the valence of the canonical transformation. So, univalent 
canonical transformations are the presymplectomorfisms between Ci and C2. 

Let {Pi, Ci, fli) be the coisotropic canonical systems associated with the 
p.l.h.s {Mi,uji,ai), i = 1,2. The class of is defined by {0} = {\E' e 
Diff(Pi,P2) I ^ ° ~ "^c o 0}- So we have the diagram 



Pi 




P2 



And therefore we have (see [7]): 



Theorem 2 There exists ip G {0} which is a symplectomorphism between 
the symplectic manifolds (Pi,f2i) and (P2,^2)- 

As a particular situation, we can analyze the canonical transformations in 
a p.l.h.s. Thus, let {M,u!,a) be a p.l.h.s., with f.c.s. jc - C ^ M. Consider 
the involutive distribution ker uc in C, and assume that the quotient space 
C = C/ker ojc is a manifold with natural projection n: C ^ C (it is called 
the reduced phase space associated to the p.l.h.s.). Then the form ujc is tt- 
projectable; hence there exists Cl G i7^(C*) such that uc = Furthermore, 
{C,Cl) is a symplectic manifold, and we have the following result (see j7]): 

Proposition 2 Every canonical transformation ($,0) in {M,uj,a) leaves 
the distribution ker uc invariant. As a consequence, there exists a unique 
G Diff(C') such that: 

1. (f) O TX = TT O (p. 

2. (j) is a symplectomorphism. 



Some applications of the geometric theory of p.l.h.s. and their canonical 
transformations are the following: the study of canonical transformations 
for regular autonomous systems (it includes a new geometrical description 
of these kinds of systems based on the coisotropic embedding theorem) [S], 
the analysis of the geometric structure and construction of canonical trans- 
formations for the free relativistic massive particle [H], the discussion about 
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time scaling transformations in Hamiltonian dynamics and their application 
to celestial mechanics jH], and the construction of realizations of symmetry 
groups for singular systems and, as an example, the Poincare realizations for 
the free relativistic particle |.35j . 

3 Constraints and the Evolution Operator 

3.1 Lagrangian Dynamical Systems 

(See PI for details). 

Let Q be a ra- dimensional differential manifold which constitutes the con- 
figuration space of a dynamical system. Its tangent and cotangent bundles, 
tq: TQ — i> Q and hq: T*Q Q, are the {phase spaces of velocities and 
momenta of the system. 

A Lagrangian dynamical system is a couple (TQ, C), where £ G C°°(TQ) 
is the Lagrangian function of the system. Using the canonical elements of 
TQ, the vertical endomorphism S G %\{TQ) and the Liouville vector field 
A G X(TQ), we can define the Lagrangian 2-form ojc, = — d(d£ o S), and 
the Lagrangian energy function Ec = A(£) — C Moreover, we define the 
Legendre transformation associated with C, T C: TQ —> T*Q, as the fiber 
derivative of the Lagrangian. 

£ is a singular Lagrangian if uc is a presymplectic form (which is as- 
sumed to have constant rank) or, what is equivalent, J-'C is no longer a 
local diffeomorphism. In particular, C is an almost-regular Lagrangian if: 
(i) Mq = J^C{TQ) is a closed submanifold of T*Q, (ii) JF£ is a submersion 
onto Mq, and (iii) for every p G TQ, the fibres jF£"^(jF£(p)) are connected 
submanifolds of TQ. Then (TQ, C) is an almost-regular Lagrangian system. 

For almost-regular Lagrangian systems, (TQ,n£,d£'£) is a p.l.h.s., and 
we have the so-called Lagrangian dynamical equation 



Variational considerations lead us to impose that, solutions Tc G X{TQ) to 
(121) must be second-order differential equations (sode); that is, holonomic 
vector fields in TQ. Geometrically this means that 



i(Tc)uJc = dEc 



(2) 



^(r^) = A 



(3) 
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For singular Lagrangians this does not hold in general, and (jH} must be im- 
posed as an additional condition (sOBE-condition) . Integral curves of vector 
fields satisfying (0) and © are solutions to the Euler-Lagrange equations. 

The Lagrangian problem consists in finding a submanifold j^^ : Sf TQ, 
and G XiTQ), tangent to Sf, such that 

[tiTc)uJc-dEc]\,^=0 , [5(r^)-A]L,^, =0 

Now, if J-'Cq : TQ Mo is the restriction of J-'C to Mq, we have that uc 
and Ec are ^£o-projectable: there exist c^o G ^^(Md), and ho G C°°(Mo) 
such that ujc = J-'CqUq, Ec = J-'CqHo. Then, (Mq, ci^o, d/io) is a p.l.h.s. which 
is called the canonical Hamiltonian system associated with the Lagrangian 
system {TQ,uJc,dEc) So we have the Hamiltonian dynamical equation 

iiXo)uJo = dho ; XoGX(Mo) 

and the Hamiltonian problem consists in finding a submanifold j^,^^ : Mf ^ 
Mo and Xq G X(Mo) tangent to Mf such that 

[i{Xo)uJo - dho]\^,^, = 
3.2 Constraint Algorithms 

In order to solve the Hamiltonian problem stated for an almost-regular sys- 
tem different kinds of Hamiltonian constraint algorithms were developed. 
The first was the local-coordinate Dirac constraint algorithm but there 
were also geometric algorithms: the Presymplectic Constraint Algorithm 
(PCA) of Gotay, Nester, Hinds [21], and others by Marmo, Tulczyjew et 
al. [Snij etc. All of them give a sequence of submanifolds which, in the 
best cases, stabilizes giving the f.c.s.: T*Q <— > Mo <— > Mi <—•...<—• Mf. 

For the Lagrangian problem, the first attempt was not to consider the 
SODE-problem (jH)), and then develop a Lagrangian constraint algorithm by 
simply applying the P.C.A. to the Lagrangian dynamical equation (j2)), ob- 
taining a sequence of submanifolds TQ <—> Pi ^ ...<—> P/ |19j . The SODE 
problem is studied later in [20]; obtaining a submanifold Sf which solves 
the lagrangian problem, but which is not defined by constraints and is not 
maximal. Later, Kamimura [2H] and Batlle, Gomis, Pons, Romdn-Roy [I] 
developed local-coordinate Lagrangian constraint algorithms in which the 
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Lagrangian dynamical equation Q and the sODE-condition (jH} were both 
considered at the same time. The f.c.s. Sf obtained at the end of the cor- 
responding sequence, T*Q S'l <—>...<—> 5/, is maximal and is defined by 
constraints. The relation between the aforementioned sequences of subman- 
ifolds is explained in the following diagram 

TQ ^ Pi ^ ... ^ Pf 

\ \ (4) 

5*1 • ... ^ Sf 

In particular, the submanifolds Pi, i = 1 . . . /, are defined by constraints 
that can be expressed as J-'C-projectable functions which give all the Hamilto- 
nian constraints, and are related with the first-class Hamiltonian constraints. 
Futhermore, the submanifolds Si, i = 1 . . . /, are defined by adding con- 
straints that are not J-'C-projectable, and they are related with the second- 
class Hamiltonian constraints. 

The remaining question was how to describe geometrically the submani- 
folds Si and their properties. First, this problem was solved in J2|; ^Sj for 
Si, i.e.; the submanifold of compatibility conditions for Q and 0, obtaining 
as the main result that: 

Theorem 3 For every SODE F G X{TQ), we have: 

Si = {peTQ\ [tiZ){t{T)uc - dEc)]ip) = 0, VZ G M} 

where M = {Z e X(Tg) | S{Z) G ker J^C, = ker ooc n X^("«)(TQ)}. 

In particular, for every Z G ker uc C A4, i{Z){i{T)uc—dEc) = i{Z){dEc) 
define the submanifold Pi where (0j is compatible. They are called dynamical 
constraints, and are related to the existence of primary first-class Hamilto- 
nian constraints. 

For every Z ^ ker Uc, Z E M., these functions are ca/Zec? SODE-constraints, 
and are related to the existence of primary second-class Hamiltonian con- 
straints. 

For Si, i = 2 . . . f , the problem was studied in j^, by imposing tangency 
conditions for solutions to (j21) and Q. The main results are the following: 
there are two kinds of Lagrangian constraints defining every Si, i = 1, . . . /: 

- Dynamical constraints, which define the submanifolds Pi, i = 1, . . . /, 
in the sequence (H)). They are related with the solutions to the eq. (0), and 
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all of them can be expressed as JF£-projectable functions which give all the 
Hamiltonian constraints. 

- SODE (non-dynamical) constraints, coming from the SODE-condition (jH}. 
They are not J^C-proj eatable. 

Hence, the relation among the sequences of submanifolds in the La- 
grangian and Hamiltonian formalisms is given in the following diagram 




T*g ^ Mo ^ Ml ^ ... ^ Mf 



Furthermore, the tangency conditions for dynamical constraints give dynami- 
cal and SODE constraints, while the tangency conditions for sODE-constraints 
remove gauge degrees of freedom and do not give new constraints. 

3.3 The Time Evolution Operator K 

As a final remark, the complete relation between Hamiltonian and Lagrangian 
constraints is given by the so-called time evolution operator K, which re- 
lates Lagrangian and Hamiltonian constraints, and also solutions to the La- 
grangian and Hamiltonian problems. (It is also known by the name of the 
relative Hamiltonian vector field It "was introduced and studied for 

the first time by J. Gomis et al. pi, developing some previous ideas of K. 
Kamimura |28j. They give the local-coordinate definition of this operator, 
whose expression in coordinates is 

The intrinsic definition and the geometric study of its properties was carried 
out independently in [10] and [23]. In the first work, K was defined using 
the Skinner-Rusk unified formalism in TQ © T*Q. In the second article, the 
concept of section along a map plays the crucial role, and so we have: 

Definition 3 Let {TQ,uJc, Ec) be a Lagrangian system, and Q G f2'^{T*Q) 
the canonical form. The time-evolution operator K associated with {TQ, uc, ^c) 
is a map K : TQ — > TT*Q verifying the following conditions: 
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1. (Structural condition) : K is vector field along T L, ttt*q o K = T L. 

2. (Dynamical condition ): T C*[i{K){Vt o C)] = dEc ■ 

3. (SODE condition): Ttq o K = Mtq- 



Then, the relation between Lagrangian and Hamiltonian constraints is 
estabhshed as follows (see [T], [TUI, ||23j): 

Proposition 3 If ^ G C°°iT*Q) is a ith- generation Hamiltonian constraint, 
then L{K)C, is a {i + l)th- generation Lagrangian constraint. 

In particular, if ^ is a first-class constraint (resp. a second-class con- 
straint^ for Mf, then L{K)C, is a dynamical constraint (resp. a SODE con- 
straint). 

The time-evolution operator has also been used for studying different 
kinds of problems concerning singular systems. For instance, the operator 
K has been extended for analyzing higher-order singular dynamical systems 
12], [H], l^n], I2Z|- It is also used for treating constrained systems in general 
[linearly singular systems) It has been defined and its properties studied 
for non-autonomous dynamical systems Finally, K has been applied for 
analyzing gauge symmetries and other structures for singular systems |22]) 
|25j . Furthermore, sections along maps in general are analyzed and used in 
different kinds of physical and geometrical problems in ^3], [T3] . 

4 Discussion and outlook 

Some of the previous problems have been studied in the sphere of first-order 
classical field theories, specially their multisymplectic formalism [3]. So, a 
geometric constraint algorithm has recently been completed for Lagrangian 
and Hamiltonian singular field theories [29^ , and the definition and properties 
of the operator K have been carried out for field theories ^7], 



TQ 



Tr, 



Q 



TT*Q 
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Other potentially interesting topics could be the generalization of some of 
the above results; such as: to study the local structure of pre-multisymplectic 
Hamiltonian field theories (previous generalization of the coisotropic embed- 
ding theorem for premultisymplectic manifolds); the study of canonical trans- 
formations for Hamiltonian field theories {multisymplectomorphisms and pre- 
multisymplectomorphisms) , and the application of the operator K to analyze 
the relation between Lagrangian and Hamiltonian constraints of singular field 
theories (which could require prior development of the non-covariant formu- 
lation, i.e., space-time splitting, of field theories). 
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